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We study the real-time dynamics of a translationally invariant quantnm spin chain, based on 
the East kinetically constrained glass model, in search for evidence of many-body localisation in 
the absence of disorder. Numerical simulations indicate a change, controlled by a coupling param¬ 
eter, from a regime of fast relaxation—corresponding to thermalisation—to a regime of very slow 
relaxation. This slowly relaxing regime is characterised by dynamical features usually associated 
with non-ergodicity and many-body localisation (MBL): memory of initial conditions, logarithmic 
growth of entanglement entropy, and non-exponential decay of time-correlators. We show that slow 
relaxation is a consequence of sensitivity to spatial fluctuations in the initial state. While nnmer- 
ics indicate that certain relaxation timescales grow markedly with size, our finite size results are 
consistent both with an MBL transition, expected to only occur in disordered systems, or with a 
pronounced quasi-MBL crossover. 


Introduction - The framework of many-body localisa¬ 
tion (MBL) [1] has been shown to be a powerful new 
approach in the study of non-equilibrium dynamics in 
closed interacting quantum systems. Since its initial pro¬ 
posal [2], abundant [3-7] numerical and analytical evi¬ 
dence has accumulated for the existence of MBL, and 
protocols for realising MBL in particular systems, such 
as ultracold atomic gases, have been proposed. MBL has 
furthermore been observed recently in an experiment on 
interacting Fermi gases [8]. While the existence of MBL 
in systems with explicit quenched disorder is now mostly 
uncontroversial, it is less clear whether a similar transi¬ 
tion can be present in disorder-free transitionally invari¬ 
ant systems [9-17]. 

Aspects of MBL have been observed in particular mod¬ 
els without disorder, such as certain mixtures of light and 
heavy particles [10, 12]. localisation in such disorder-free 
systems would occur due to the interplay between inter¬ 
actions and inhomogeneous initial conditions, leading to 
non-ergodic behaviour where relaxation times would di¬ 
verge with system size. In contrast, a general mechanism 
for delocalisation has been suggested in [15], by which 
the presence of mobile resonant spots act as carriers of 
energy or charge. At present it is unclear whether true 
MBL can occur in a finite-dimensional translation invari¬ 
ant system, and quasi-MSL has been suggested [12] for 
cases where transient MBL behaviour is followed by the 
restoration of ergodicity on long timescales. (The situa¬ 
tion is reminiscent in some sense to that of the classical 
glass transition problem, where it is debated whether re¬ 
laxation times diverge at finite temperature with a tran¬ 
sition to an ideal glass state, or whether relaxation times 
are very long but remain finite for all non-zero tempera¬ 
ture. See e.g. the discussion in Ref. [18].) 

In this work we argue for a transition between an er- 
godic (delocalised) phase and an MBL phase in a specific 
translation invariant system. Generally, a disorder less 
delocalised-to-MBL phase transition is conjectured to be 
related to a quantum glass transition [1]. In particu¬ 
lar, in Ref. [11] indications of such a delocalised-to-MBL 


transition were found in a one-dimensional quantum spin 
chain based on the classical glassy Fredrickson-Anderson 
model [19, 20]. Here we study an even simpler system— 
and for which MBL-like dynamics is even more marked— 
a disorderless quantum spin chain based on the (infinite 
temperature) Fast glass model [20-22]. We find that in 
this model a single coupling parameter controls a sig¬ 
nificant change in the dynamics, from a regime of fast 
relaxation and thermalisation, to one where some relax¬ 
ation times appear to diverge exponentially with system 
size. We show that in this slowly relaxing regime the sys¬ 
tem exhibits many of the dynamical features associated 
with MBL: logarithmic spreading of excitations, lack of 
thermalisation, initial state dependence at long times and 
logarithmic growth of entanglement entropy. While these 
observations are highly suggestive of MBL we cannot dis¬ 
card quasi-MBL behaviour as an alternative explanation. 

Model- We consider a translation invariant spin-1/2 
chain with Hamiltonian operator 

N N 

H = -^aluk+i. ( 1 ) 

fc=l k=l 

Here J is an energy scale, which, without loss of general¬ 
ity, we take to be J = 1. On the other hand, varying the 
parameter s, which controls the coupling between sites, 
drastically changes the dynamical behaviour of the model 
as can be observed in Fig. 1. This Hamiltonian is closely 
related to the master operator of the classical East glass 
model [20-22]. In fact, if used to generate imaginary time 
dynamics, the s = 0 Hamiltonian in Eq. (1) is equivalent 
to the master operator of the classical East model at infi¬ 
nite temperature [i.e., s = 0 is a Rokhsar-Kivelson point 
[23, 24] of (1)]. At s 7 ^ 0 the operator H is a so-called 
tilted generator [25], whose s-dependent ground state is 
the cumulant generating function of the “dynamical ac¬ 
tivity” [26]. The ground state of the Hamiltonian (1) is 
known to change discontinuously at s = 0 [27]. (See also 
the discussion in [11] for a related model.) 

Classical mechanism for relaxation - In the classical 
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FIG. 1. (Colour online) Time evolution of site densities {ni{t)} 
starting in the half-filled initial state |4, • • • • • • t) in 

the thermalising regime, s = —0.5 (a), and in the MBL-like 
regime, s = 0.6 (b), showing a linear (ballistic) spreading of 
the state and a logarithmic spread, respectively, (c): equilib¬ 
rium value Xeq of the initial participation ratio T{t), averaged 
over all half-hlling initial states. 


stochastic dynamics generated by (1) when s = 0, con¬ 
figurations are represented by Fock states in the 2 ba¬ 
sis. The key feature of these dynamics is that changes 
of configuration occur via individual spin-flips only on 
sites which have a right nearest neighbour with cr^ = 1. 
This means that a region of space with contiguous down 
spins can only be relaxed from its rightmost boundary, 
and as a result, longer regions of this type take longer to 
relax (as the connecting dynamical pathway is required 
to show more collective features). This “rare region” be¬ 
haviour is what makes the East model glassy at low tem¬ 
peratures (when the “facilitating” up spins are sparse) 
[20-22]. Furthermore, fluctuations in the initial config¬ 
uration associated to these down spin domains lead to 
fluctuations in the spatial relaxation pattern—known as 
“dynamical heterogeneity”—of the classical East model 
[28]. As we will see below, analogous initial state fluc¬ 
tuations strongly affect the coherent dynamics generated 
by Eq. (1) for general s. 

Dynamical features - Considering pure initial states 
which lack translational invariance reveals dynamical fea¬ 
tures of our model when unitarily evolving the system 
with the Hamiltonian H. As the system is closed and 
subject to recurrence one has to consider time-integrated 
observables to account for relaxation, or limit the dynam¬ 
ics to a small portion of the recurrence cycle [29, 30]. 

We find two distinct dynamical regimes parametrized 
by s. For s < 0 there is fast relaxation and ballistic 
spreading of initial excitations (up spins), as illustrated in 
Fig. 1(a). In contrast, the s > 0 regime exhibits slow re¬ 
laxation and logarithmic spreading of initial excitations, 
as shown in Fig. 1(b). In both figures the initial state is 
a half-filling Fock state (a Fock vector with equal num¬ 
ber of up and down spins in the z basis) where the first 
half of spins are down and the rest is up. We now argue 
that the two regimes are distinguished by thermalisation 



FIG. 2. (Colour online) (a,b): time-averaged magnetization 
Mn,{t) resolved by the largest interval w between excitations 
in the initial state (c): initial state dependence quantified by 
var[M]eq/var[M]niax (d): comparison between the equilibrium 
magnetization Meq, and its microcanonical average (M)mc. 
The highlighted points are the equilibrium values Meq ob¬ 
tained in panel (b), corresponding to initial states with energy 
E « 1.7. 

when s < 0 and lack of thermalisation for s > 0, and 
the non-thermalising regime in particular exhibits many 
of the features of a MBL system. 

To probe the dynamics for the possible occurrence of 
a thermalising-to-MBL transition we consider the dy¬ 
namical inverse participation ratio (IPR) [31] I{t) = 
2~^ (Si l(V'(^)|i) I"*) ^ associated to the state at 

time t > 0. Fig. 1(c) shows the equilibrium [32] value of 
the IPR leq, averaged over all half-filling initial states. 
The quantity leq acts as an indicator of localisation of 
the equilibrium state in the Fock basis; in particular, 
its value changes from approximately unity for s < 0 to 
0{N~^) for s > 0, suggesting that at the point s « 0 a 
MBL transition might occur. 

Initial state dependence - We consider the dynamical 
features of the system starting from a pure initial state 
picked randomly from the half-filling sector. We charac¬ 
terise the initial state dependence by associating to each 
initial classical spin configuration a quantity w, defined 
as the maximal number of consecutive down spins, taking 
into account periodic boundary conditions (e.g. w = 3 
for a Fock vector jitiittitii))- 

Long-time dependence of single realizations on the ini¬ 
tial state indicates a lack of thermalisation, and is con¬ 
sidered a candidate indicator of many-body localisation 
[3, 4, 33, 34]. We find that the equilibrium state in the 
regime s > 0 shows non-thermalising behavior by retain¬ 
ing memory of initial conditions even at long timescales. 
This can be seen by considering the time-integrated total 
magnetization in the z-direction M{t) = t~^ m{t')dt', 
where m{t) = (<p(<)| Fig. 2(a-b) we dis¬ 

play this quantity respectively in the thermalising and 
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FIG. 3. (Colour online) (a): entanglement entropy S{t) (with respect to an equal bipartition of the spin chain) showing 
qualitative difference as s is varied. Inset: equilibrium value S'eq of entanglement entropy as function of s and N, suggesting 
a transition from the thermalising regime s < 0 to the MBL-like regime s > 0. (c,d): average entanglement entropies Sw{t) 
restricted to initial states with domain size w, showing that initial state dependence is absent in the thermalising phase, (d): 
Characterization of the logarithmic regime in the growth of the entanglement entropy S{t). Dashed lines shows best logarithmic 
fit S{t) oc logt for each value of N. Inset: estimated relaxation time r associated to S{t). 


non-ergodic regimes resolved according to the maximum 
distance of excitations present in the initial state Fock 
vector. We denote with Mw{t) the average over realiza¬ 
tions restricted to initial states with maximum distance 
of w sites between excitations. 

In the regime s < 0 the system does not retain 
memory of the initial state in the long-time limit, as 
Mw{t) 0 for all values of w. When s > 0 on 
the other hand the long-time values of Mw{t) show 
a clear dependence on w, suggesting an equilibration 
which depends on the initial conditions, and as such 
no thermalisation occurs. This lack of thermalisation 
is reflected in the time-dependent variance (with re¬ 
spect to initial state realizations) var[M(t)]. In Fig 2(c) 
we make the potential transition apparent by plotting 
the quantity var[M]eq/var[M]i„ax 7 where var[M]i„ax = 
maxoo var[M(t)] and var[M]eq denotes the equilibrium 
value of var[M(t)]. This scaled variance compares the 
maximum transient spread of the integrated magnetiza¬ 
tion (a variation which is expected since at intermediate 
times different initial states will evolve differently) to the 
spread at long times (a variation which only occurs in 
the absence of thermalisation). 

Further evidence for the absence of thermalisation for 
s > 1 is shown in Fig. 2(d) where we plot the equilib¬ 
rium values of the magnetization Mgq along with its mi- 
crocanonical average expectation value (at infinite tem¬ 
perature), as a function of energy [35]. Here it 

can be seen that the microcanonical prediction does not 
reproduce the long-time average when the initial state 
energy lies in the bulk of the spectrum. In this case Meq 
strongly depends on the initial conditions; these results 
suggest a breakdown of the eigenstate thermalisation hy¬ 
pothesis leading to absence of thermalisation [36], which 
in turn is connected to MBL [3]. 


Logarithmic growth of entanglement entropy- Slow 
growth of the entanglement entropy has been associated 
to a MBL phase in both translation invariant [12] and 
disordered systems [4, 37]. We consider the entangle¬ 
ment entropy with respect to an equal splitting of the 
system into two half chains, denoted A and B. Let 
S{t) = —Tr^pyi(t) log pA(i) be the time-dependent en¬ 
tanglement entropy, where pA{t) = ATB\^{t)){if{t)\ is 
the reduced density matrix on the half chain A. We de¬ 
fine S(t) as the mean entanglement entropy obtained by 
taking the average over the S{t) associated to all pos¬ 
sible initial states in the half Filing sector; additionally, 
to maintain translation invariance, S(t) is also averaged 
over all possible locations of the cut between half chains 
A and B. 

In Fig. 3(a) we compare the evolution of the average 
entropy S{t) for different values of s, starting from the 
same initial half-filling ensemble. Clearly evident is a dif¬ 
ference of behavior in the thermalising (s < 0) and the 
non-thermalising (s > 0) regimes: the former is char¬ 
acterized by an exponential growth which culminates in 
the equilibrium entropy S'eq which is independent of s. In 
the non-thermalising regime on the other hand, after an 
initial fast growth we observe a slow logarithmic increase 
which extends over several order of magnitude of Jt. The 
equilibrium value Seq in this regime drops significantly 
for increasing s, as shown in the inset in Fig. 3(a). 

The time evolution of S(t) can be understood quali¬ 
tatively as follows: As previously discussed the typical 
initial state will contain a largest domain of down spins 
of length w. The initial linear increase of S{t) corre¬ 
sponds to the mixing of the state in the regions in away 
from this empty interval. Once the state is completely 
mixed over these regions the only increase in entangle¬ 
ment can come from propagating excitations into the w 
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down spin domain, which as shown in Figs. l(a,b) is lin¬ 
ear in time for s < 0 and logarithmic for s > 0. In the 
thermalising regime, s < 0, there is no difference in the 
growth of the entanglement entropy starting from states 
with different w, and equilibration is fast. In contrast, in 
the non-thermalising regime, s > 0, this growth shows a 
pronounced dependence on the value of w in the initial 
state. 

The range of logarithmic growth of the entropy in the 
s > 0 regime increases with system size, as shown in Fig. 
3(d). We can extract a relaxation time r as the time the 
S{t) reaches Seq, for example such that | Seq-S{T)\=e 
with e = 0.01 (the precise value of this cutoff is unimpor¬ 
tant) [38]. The inset to Fig. 3(d) we shows the scaling of 
the relaxation time r for different values of s. For s > 0, 
T it appears to grow exponentially with the system size 
in the range of sizes which we were able to consider. This 
exponential divergence with size is a hallmark of the loss 
of ergodicity and of MBL. 

Polarization relaxation — Finally, we consider the re¬ 
laxation of fluctuations at different wavelengths. To do 
this we consider the (infinite temperature) averaged time- 
correlator of the Fourier transformed longitudinal mag¬ 
netization, D{k,t) = Tr[e“*^‘F*(fc)e*^‘F(fc)], where 
F{k) = The relaxation of this polarization 

correlator can be used to quantify the transport of spin 
perturbations [3, 12] and thus can probe for MBL. Since 
in a translation invariant system of finite size, D{k,t) 
has to decay to zero at large enough times, the interest¬ 
ing property to look for is in the (transient) behaviour of 
the relaxation times [12]. 

In Fig. we show D{k, t) as a function of time for all the 
wave vectors fc of a system of size TV = 12 (the largest 
we were able to simulate for this kind of correlator). For 
the thermalising case of s = —1, Fig. 4(a), all correla¬ 
tors decay exponentially and there is no fc-dependence 
in their relaxation. In contrast, for the non-thermalising 
case of s = 1, Fig. 4(b), relaxation is non-exponential and 
timescales are fc-dependent. Figure 4(c) shows estimates 
of the relaxation times t as a function of k for three sys¬ 
tem sizes: for s < 0, r(fc) is independent of k (and of 
system size), while for s > 0, they are strongly k depen¬ 
dent. In fact, in the non-thermalising phase T{k) ~ 
for a large range of wavelengths, a similar behaviour as 
that observed in Ref. [12]. 

The exponential in wavelength behaviour of the relax¬ 
ation time can be rationalized as follows. D{k,t) mea¬ 
sures the relaxation over a wavelength 2TTjk. In the MBL- 
like regime, the initial states that couple to this wave¬ 
length and decay the slowest [and thus dominate D{k,t)] 
will be those with w = 2TT/k, where w is the size of the 
largest domain of down spins. As shown above, such do¬ 
mains decay logarithmically in time, cf. Fig. 1(b), that is 
the relaxation into an inactive domain goes as ru « ^ log t. 
Inverting this relation and replacing w hy 1/k we get the 
form of T(fc) above. In order for this argument to hold 



FIG. 4. (Colour online) Dynamics of the spin polarization 
D{k, t) for all wave vectors fe in a system of size N = 12. 
(a): D{k,t) in the thermalising phase, s = —1. (b): D{k,t) 
in the non-ergodic phase, s = 1. (c): Wavelength dependent 
relaxation timescales, estimated by D{k,T) = e~^, at three 
values of s; the dashed lines indicate scaling. 

these slow initial states have to be numerous in the infi¬ 
nite temperature average. At the largest wavelength they 
become rare, and this is be the reason why the exponen¬ 
tial scaling is not obeyed in Fig. 4(c) for small k. 

Conclusions and outlook - In this paper we studied a 
disorderless quantum spin chain closely related to the 
classical glassy East model to explore a possible MBL 
transition in a translationally invariant system. The co¬ 
herent dynamics of this model displays two regimes; a 
thermalising regime, and a non-ergodic regime charac¬ 
terised by strong dependence on the initial conditions 
and breaking of thermalisation. The non-thermalising 
regime displays dynamical features characteristic of an 
MBL phase, in particular a logarithmic spreading of ex¬ 
citations and a logarithmic growth of the entanglement 
entropy. The slow relaxation behind these MBL like fea¬ 
tures can be traced back to spatial fluctuations in the ini¬ 
tial state corresponding to rare inactive domains. While 
relaxation times seem to diverge with system size for the 
entanglement entropy, they remain finite (yet anomalous) 
for dynamic structure factors. These features could be 
compatible both with MBL [1] or with quasi-MBL [12]. 
In either case, our results suggest a close connection be¬ 
tween classical mechanisms for glassy arrest and those 
responsible for MBL dynamics. 
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